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We have investigated the magnetic phase diagram of a fully oxygenated detwinned YBa2Cu307 
single crystal by means of magneto-caloric and magnetization measurements, and found thermo- 
dynamic evidence for a temperature dependent first-order phase-transition line deep within the 
vortex-solid region. The associated discontinuities in the entropy are apparently proportional to the 
magnetic flux density, which may hint at a structural transition of the vortex lattice. 
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The vortex lattice of type-II superconduetors in the 
mixed state can undergo phase transitions. Besides the 
well known melting transition in cuprate superconduc- 
tors [1-10], other possible kinds of transitions of the vor- 
tex lattice have become of growing interest. Nonlocal 
effects in the vortex-vortex interaction that are present 
in certain low-temperature superconductors such as the 
borocarbides or VsSi lead to a rich magnetic phase dia- 
gram, with, e.g., structural phase transitions from cu- 
bic to hexagonal vortex lattices [11, 12]. The struc- 
ture of the vortex lattice and its phase transitions have 
been extensively investigated in the past using differ- 
ent imaging techniques [13-27]. However, thermody- 
namic data on such structural vortex phase transitions 
(other than the melting transition) are extremely rare 
[28]. The reported first-order transitions in the specific 
heats of CeCoIns [29] and in the organic superconductor 
K-(BEDT-TTF)2Cu(NCS)2 [30] have been interpreted 
in terms of a Fulde-Ferrell-Larkin-Ovchinnikov scenario 
that invokes a vortex segmentation, rather than a mere 
geometrical rearrangement of the vortex lattice. 

Very recent small-angle neutron-scattering (SANS) 
experiments on YBa2Cu307 have revealed Fermi sur- 
face and/or order-parameter driven, previously unknown 
structural phase transitions deep within the vortex-solid 
region, i.e., well below the upper-critical field Hc2{T) 
and the vortex-lattice melting field H,n{T) [26]. At 
the temperature T = 2 K and with the magnetic field 
H applied along the c-axis, the vortex-lattice structure 
changes abruptly at (J-qH w 2T from a low- field hexag- 
onal to an intermediate-field distorted hexagonal phase, 
and at iiqH « 6.7 T to a high-field rhombic structure. 
On the basis of the observation of discontinuous changes 
in the lattice structure, these transitions have been sug- 
gested to be of first order and the transition line was 
found to be temperature independent [26]. 

In this work we report on the observation of a tem- 
perature dependent first-order phase transition line deep 
within the vortex-solid region. We conducted measure- 
ments of the quasi-isothermal magneto-caloric effect in 



an ultrapure YBa2Cu307 single crystal [31] that is part 
of the mosaic studied in these SANS experiments. The 
crystal has a critical temperature Tc ~ 86 K and in a 
magnetic field of 7T the melting of the vortex lattice 
was observed at T™(7T) w 80K [32]. To reduce vortex 
pinning on defects, the crystals had been slightly over- 
doped to achieve full stoichiometry, and they were also 
mechanically detwinned [26, 33]. The crystal used in our 
experiments (mass ~ 6.5 mg) was aligned with its c-axis 
parallel to the external magnetic field. The calorimet- 
ric measurements were done in a home-built differential- 
thermal analysis (DTA) setup [34] and run in a constant- 
temperature mode with varying magnetic field. In exper- 
iments of this type, the sample and its suspension with 
total heat capacity C are connected to a thermal bath 
with a known heat link k that defines the time constant 
T = C/k according to which the sample can approach 
the temperature of the thermal bath. On varying the 
external magnetic field from a starting value Hq to H 
at a rate dH/dt, a deviation AT{H) of the sample tem- 
perature from the constant thermal bath temperature T 
evolves that can be used to calculate the variation in the 
entropy S{H) according to 

[S{H) SiHo)]^ « - [AT{H) - AT(ffo)] (1) 



Jh, r dH/dt 

which is valid for AT -C T, i.e., for small variations 
of S as compared to C [34]. While the first term 
in equation (1) originates from the definition of the 
entropy S, the second one represents a deconvolution 
term accounting for the effect of a finite heat link on 
the sample temperature in such an experiment, and it 
vanishes under perfectly adiabatic conditions, i.e., fc = 0. 

White et al. [26] prepared the vortex lattice for the 
above cited SANS measurements according to two dif- 
ferent techniques. The first one was a standard field 
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cooling through Tc in a constant magnetic field ('field 
cooling', FC). The second one was a field cooling with 
the main (longitudinal) field oscillating around the tar- 
get value ('oscillation field cooling', OFC) with an os- 
cillation amplitude of the order of a few mT, and at 
a frequency / « 30mHz. For both FC and OFC, the 
field was held stationary for diffraction measurements at 
T = 2K. The resulting SANS patterns obtained from 
OFC measurements showed a significant improvement 
of the vortex-lattice perfection as compared to measure- 
ments conducted without using an oscillating field. 

In most of our measurements we cooled the sample 
to the target temperature in a constant magnetic field 
/io-ff = 7T prior to a measurement. The magneto-caloric 
effect was then monitored while decreasing H to zero at 
a constant rate ^odH/dt w — 1 mTs^^. Since our mag- 
net power supply is not able to ramp the magnetic field 
perfectly linear, an ac magnetic field component paral- 
lel to the main magnetic field was always superimposed 
in such experiments with an amplitude /^o^ac 

« 0.4 mT 

and a frequency / « 11 mHz, resulting in conditions very 
similar to those for the OFC technique used in [26]. Ac- 
cording to [35] the magnetic field of full penetration is 
Bp ~ fiojcW, where 2w is the width of the sample and 
jc the critical-current density. With w « 0.75 mm and 
jc « 2.88 X 10^ Am^2 at T = 10 K that we estimated from 
our magnetization data, we obtain Bp « 2.7 T, which is 
much larger than the oscillation amplitudes used in both 
the SANS study and in our experiments, so no full pen- 
etration of the ac magnetic field can be expected. 

In figure 1(a) we show the resulting magneto-caloric 
effect for various temperatures on ramping the magnetic 
field H II c down (i.e., with Hq > H and dH/dt < in 
equation (1)). (The (5 AT data represent the AT data 
after subtracting a polynomial background.) Abrupt in- 
creases in AT{H) at certain magnetic fields H*{T) indi- 
cate the occurrence of a first-order phase transition sep- 
arating a higher-field phase from a low-field phase with 
lower entropy [34]. These AT{H) data can be deconvo- 
luted according to equation (1), and an example of the 
resulting background corrected entropy change is shown 
in figure 1(b). In figure 2(a) we have plotted the cor- 
responding magnetic phase diagram, together with the 
extracted discontinuities in the entropy AS" in an inset 
to that figure. The experiments below 10 K were too 
noisy for extracting reliable AS values due to unwanted 
thermal effects induced by the superimposed ac magnetic 
field, although the transition at H* was still clearly dis- 
cernible also in these AT{H) data. 

If the reported features were of thermodynamic ori- 
gin they ought to be seen also in experiments in the op- 
posite direction, i.e., when crossing the supposed phase 
boundary H*{T) from the low-field to the higher-field 
phase. In related experiments in a constant magnetic 
field (i.e., without any oscillating field component) but 
with increasing temperature, we did indeed observe a cor- 
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FIG. 1: (a) Magneto-caloric effect measured in an 
YBa2Cu307 single crystal on ramping the magnetic field 
II c down. A polynomial background lias been subtracted 
to obtain SAT, and the curves are vertically shifted for clarity. 
The large red circles mark the field which we defined as the 
transition field H* . (b) Example illustrating the procedure 
used to extract the changes in the entropy according to [34] . 



responding cooling effect on crossing the phase-transition 
line H*{T) from below (see figure 2(b)). A correspond- 
ing data point is also drawn in the phase diagram in fig- 
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FIG. 2: (a) Phase diagram of a detwinned stoichiometric 
YBa2Cu307 single crystal as derived from magneto-caloric 
effect and magnetization data. The SANS data point from 
[26] has been added to the phase diagram for comparison. 
The line is an exponential fit to the magneto-caloric data 
(see the text). The grey and orange shaded areas are derived 
from the transition widths of the magneto-caloric and magne- 
tization data, respectively. Inset: entropy difference between 
the higher-field and the low-field phase. The line is a linear 
fit. (b) DTA temperature scan in a constant magnetic field 
fioH = 0.7 T after subtraction of a polynomial background. 
The data are showing the cooling effect due to the latent heat 
absorbed at the first-order phase transition. 



ure 2(a). Note that the magnitude of SAT in the latter 
type of measurement (figure 2(b)) is smaller than that 
measured in the field-sweep experiments shown in fig- 
ure 1 , which is a direct consequence of the relatively low 
measuring speed with which the phase-transition line has 
been crossed. According to [34], corresponding signals in- 
volving the same latent heat TAS scale with the inverse 
of the time needed to pass a phase transition of finite 
width (see below) , which was w 750 s for the temperature 
scan and « 72 s for the field scan at approximately the 
same temperature T « 16.5 K. The observation of both 
heating and cooling effects, depending on the direction in 
which H* (T) is crossed, strongly supports our interpreta- 
tion that the measured variations in sample temperature 
are in fact related to abrupt changes in the entropy of 
the magnetic sub-system (i.e., to the release or the ab- 
sorption of a latent heat to or from the crystal lattice, 
respectively). They can therefore not simply be caused 
by dissipative thermal effects, e.g. induced by a possible 
vortex motion upon the rearrangement of the vortex lat- 
tice, that would lead to a heating effect in all cases. 

In figure 3(a) wc show magnetization data M{H) 
measured using a SQUID magnetometer (Quantum De- 
sign Inc.) for the same single crystal as we have used 
in the thermal experiments. On ramping the magnetic 
field down, we observed step-like features in M that are 
reminiscent of a first-order discontinuity AM and oc- 
cur at magnetic fields which essentially coincide with the 
phase boundary H*{T) as obtained from the magneto- 



caloric experiments (see figure 2(a)). However, these 
features are embedded in a background magnetization 
with gradual and perhaps even discontinuous changes 
in the slope dM/dH around the proposed phase tran- 
sition at H* . This fact and the presence of a size- 
able irreversible contribution to M (see figure 3) ren- 
der a reliable extraction of a AM somewhat difficult 
as the result may depend on the assumptions to de- 
scribe the underlying background magnetization. In an 
attempt to make the background-suppression procedure 
as unbiased as possible, we have calculated the aver- 
aged magnetization Mav{H) from M{H) data taken for 
decreasing and increasing magnetic field, respectively, 
and we plotted the resulting values in the region of 
interest in figure 4. For a tentative comparison with 
the magneto-caloric data, we have also drawn in esti- 
mates for a hypothetical equilibrium discontinuity AM 
that we obtained from the Clausius-Clapeyron equation, 
AS{H*) = -AMnodH*/dT, where we used AS" values 
from a linear fit to the thermal AS{H*) data shown in the 
inset of figure 2(a). Although the magnetization is hys- 
teretic and corresponding step-like features arc virtually 
absent in the M{H) data taken for increasing magnetic 
field (see figure 3(b)), the agreement of the AM values 
estimated from the magneto-caloric experiments with the 
order of magnitude of the measured step-like features in 
Mav, as shown in figure 4, is quite fair, and it may hint 
at the occurrence of a true, thermodynamic phase tran- 
sition at H*{T). 

For completeness we would like to mention that the 
field inhomogeneity of the solenoid used, providing the 
main magnetic field over the used scan length (4 cm) used 
in the DC mode leads to an unavoidable superimposed 
longitudinal ac magnetic field of the order of fiohac ^ 1- 
3 mT at a leading measuring frequency of ~ 30 niHz that 
must have been present during all of these magnetization 
measurements. 

The width of the phase transition at T « 17K is 
fioH*^^ « 70 mT as derived from the magneto-caloric field 
scans shown in figure 1(b). This is comparable to the re- 
gion where the step-like features and the slope changes in 
the M (H) data occur (see the horizontal arrow in figure 
4). It also coincides with the width T^^ « 1.5 K obtained 
from the temperature scan shown in figure 2(b) that can 
be converted to /ioHt*^ ~ T^^^^QdH* /dT w 75 mT using 
the slope fit^dH* /dT = — 0.05TK~^ at this temperature. 
This good agreement of the transition widths lends sup- 
port to our assumption that the transition observed in 
the temperature sweep is the same transition as was ob- 
served in the field sweeps. 

The phase-transition line, as constructed from 
our thermodynamic measurements, fairly accurately 
follows the phenomenological relation fiQH*{T) = 
2.3 exp(— r/13.5 K) T and extrapolates at low tempera- 
tures to a data point where detailed SANS experiments 
taken at T = 2 K and in /loi? ~ 2 T have revealed a trans- 
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FIG. 3: Magnetization M{H) of the investigated YBa2Cu307 
single crystal for (a) decreasing and (b) increasing magnetic 
field II c for different temperatures. To make the region of 
the proposed phase transition at H* clearer (see the arrows), 
we have shifted the curves horizontally and vertically on the 
same relative scales, and indicated the absolute values of fioH 
and M for each M{H) curve at a selected data point. The 
straight lines are fits to the AI (H) data for H < H* . 




FIG. 4: Averaged magnetization Mav{H) (see the text), 
shifted in an analogous way to that in figure 3. Next to the 
arrows which mark the phase transition at //* we plotted ver- 
tical bars representing estimates for the expected step in the 
magnetization as derived from our magneto-caloric data and 
the Clausius-Clapeyron equation (see the text). The horizon- 
tal arrow represents the transition width fioHt^ as determined 
from the magneto-caloric measurements (see figure 1(b)). 



formation of the vortex lattice from a low-field hexago- 
nal to an intermediate-field distorted hexagonal phase 
[26] (see figure 2(a)). However, additional SANS mea- 
surements of the vortex lattice form factor at five dif- 
ferent temperatures point to a virtually constant phase- 
transition field up to the vortex-lattice melting line [36]. 
This obvious discrepancy to the phase diagram shown 
in figure 2(a) may indicate that the phase transition de- 
scribed here is different from that reported in [26] and 
[36], and it may therefore also have a different physical 
origin. The phase transition is probably not a kinetic 
glass transition from an ordered to a disordered vortex 
glass since such a transition was recently reported by 
Petrovic et al. [27] to be not accompanied by a latent 
heat. A relation to an order-disorder transition in the 
vortex solid as observed in [37] is also unlikely since the 
phase-transition line as derived from our experiments is 
located at much lower temperatures and magnetic fields 
in the phase diagram. 

We want to emphasize that the magneto-caloric mea- 
surements probe a bulk property of the ensemble of vor- 
tices whose number is given by the flux density B, and 
we therefore judge our magneto-caloric AS' data as more 
reliable than corresponding estimates from the magne- 
tization HqM <^ _B, which is highly susceptible to irre- 
versible behavior. It is noteworthy that the measured 
AS{H*) scales apparently linearly with the vortex den- 
sity (see the inset of figure 2(a)). However, we obtain a 
AS ~ 20 fcs per vortex per CU-O2 double layer (with the 
Boltzmann constant ks), which is fairly large for a purely 
structural phase transition of the vortex lattice where one 
might expect AS" to be of the order of 1 fc^ per vortex per 
double layer. Therefore wc do not rule out, as already 
mentioned, the possibility that the observed first-order 
phase transition may have an entirely different physical 
origin, and that the apparent agreement of the extrap- 
olated phase-transition line H*{T) with the SANS data 
at T = 2 K is coincidental. 

In summary, we conclude that we have found a thermo- 
dynamic evidence for a true first-order phase transition 
deep within the vortex-solid region of YBa2Cu307. The 
associated discontinuities in the entropy scale virtually 
linearly with the magnetic field, which may indicate that 
we observe a structural phase transition of the vortex 
lattice. Since the transition is of first order, a glass tran- 
sition can be ruled out as a possible origin. However, the 
entropy changes are an order of magnitude larger than 
what has been measured at the known vortex solid-to- 
liquid transitions, and we can therefore not rule out an 
entirely different physical origin of the observed phase 
transition. 
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